High-energy photoemission on Fe 3 04: Small polaron physics and the Verwey 

transition 
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We have studied the electronic structure and charge ordering (Verwey) transition of magnetite 
| (FesCU) by soft x-ray photoemission. Due to the enhanced probing depth and the use of differ- 

$_h , ent surface preparations we are able to distinguish surface and volume effects in the spectra. The 

pseudogap behavior of the intrinsic spectra and its temperature dependence give evidence for the 
existence of strongly bound small polarons consistent with both dc and optical conductivity. To- 
gether with other recent structural and theoretical results our findings support a picture in which the 
Verwey transition contains elements of a cooperative Jahn- Teller effect, stabilized by local Coulomb 
interaction. 



PACS numbers: 71.30.+h, 71.38.Ht, 79.60.-i 



Magnetite (Fe3 04) is not only the oldest magnetic material known to mankind, with high potential 
for applications in spin-electronics, but also displays a rather unique electronic phase transition whose 
explanation has remained a challenge to modern condensed matter physics 0, 0] • At low temperatures 
. magnetite is an insulator. Upon heating its resistivity drops abruptly in a first-order phase transition 
at Ty ~ 123 K by about two orders of magnitude to that of a bad metal, although its temperature 
. dependence up to 300 K remains non-metallic. Verwey Q was the first to suggest that the transition 
is related to charge ordering (CO) in this mixed-valent oxide. In this picture the Fe 3+ and Fe 2+ ions, 
which occupy the octahedral B-sites of the inverse spinel structure in equal numbers, form below Ty an 
ordered pattern, whereas above the transition they become statistically distributed over the _B-sublattice 
(with the possible preservation of short-range order). The driving force for the Verwey transition has 
been suggested to be the interatomic Coulomb interaction [4|. However, as first pointed out by Ihle and 
Lorenz strong coupling to the lattice may also be of importance. 

While the Verwey transition has thus been studied for decades, very recent experiments have raised 
fundamental questions. Based on new high-resolution powder diffraction data it has been estimated that 
the actual charge disproportionation is only very small, of order ±0.1e 0|. From resonant x-ray diffraction 
other authors even concluded on the complete absence of any CO in the insulating phase of magnetite 
0, @- These results suggest that the electrostatic energy minimization implied by the proposed CO 
cannot be the dominant force behind the Verwey transition contrary to previous ideas. 

Another controversial question concerns the electronic nature of the high-temperature phase. In a high- 
resolution photoemission study Chainani et at [jj found above Ty finite spectral weight at the chemical 
potential and concluded on a metallic state. In contrast, Park et al. [13 reported photoemission spectra 
measured just below and above Ty which display a notable change in the gap but leave it finite above 
Ty. The Verwey transition was thus identified as insulator-to- insulator transition. We note however that 
all previous photoemission experiments have been performed at low photon energies, where the probing 
depth 11] is small (~ 15 A) and the spectra are very surface-sensitive. 

In this letter we demonstrate that low-energy photoemission spectra of magnetite are indeed strongly 
affected by surface orientation and preparation, thus limiting their reliability for a study of the bulk 
electronic structure. We show how this problem can be overcome by use of soft x-ray photoemission 
which due to its enhanced probing depth (~ 45 A) allows much better access to the volume. It is 
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observed that the intrinsic spectra both above and below Ty show an exponential suppression of spectral 
weight towards the chemical potential and bear no resemblance to the prominent doublet of quasiparticle 
peak and lower Hubbard band as observed in other oxides with Mott-like metal-insulator transitions. In 
fact, from the temperature dependence of the spectra we conclude on a semiconductor-semiconductor 
character of the Verwey transition. We analyze the pseudogap-like line shape and demonstrate that 
it can consistently be attributed to strong electron-phonon coupling. The results yield evidence that 
small polaron physics plays an essential role for the electronic properties and the Verwey transition of 
magnetite. 

The experiments were performed on high-quality synthetic magnetite crystals. From these we cut ori- 
ented (111) and (100) surfaces which after polishing were exposed to an in situ treatment involving Ar ion 
sputtering and annealing (up to 800° C). Loss of surface oxygen was compensated by a subsequent oxida- 
tion step (further details are described in 0). Long-range order and stoichiometry of the surfaces were 
checked by LEED and XPS, respectively. Alternatively, clean surfaces were also obtained by fracturing 
single crystals in situ at 100 K. Because magnetite has no natural cleavage plane, the resulting surfaces 
are rather rough and faceted. Low-energy photoemission was performed at our home lab u sing He I 
radiation (hv — 21.2 eV). The soft x-ray experiments were performed at BL25SU of SPring-8 [13j using 
Fe 2p-3d resonance photoemission (hv — 707.6 eV). In both cases the energy resolution was 100 meV or 
lower. The position of the chemical potential was calibrated by the Fermi edge of an evaporated Au film 
in electrical contact with the magnetite crystals. 
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FIG. 1: (a) Valence band spectra of three differently prepared and oriented magnetite surfaces, measured with 
vacuum ultraviolet radiation. Inset: Spectral weight distribution near the chemical potential. Note the dramatic 
surface effects on the low-energy electronic structure, (b) Corresponding spectra measured with soft x-ray ra- 
diation at the Fe L2 edge. The strongly reduced differences between the various surfaces give evidence of the 
enhanced bulk sensitivity. 

Figure f^a) an d (b) show the valence band spectra of the different magnetite surfaces measured with 
vacuum ultraviolet (VUV) and soft x-ray (SX) radiation, respectively. The VUV spectra display dra- 
matic differences depending on surface preparation and orientation, reflecting strong electronic structure 
rearrangements in the outermost atomic layers. In contrast, the SX spectra show only very small (though 
reproducible) differences, indicating that due to the larger probing depth they are dominated by signal 
from the bulk. This is particularly evident in the spectral weight distribution near the chemical potential 
/x (see insets in Fig. ^) which originates from the d 6 — > d 5 transition at the B-site Fc 2+ ions 0, 0] and is 
directly related to the low-energy physics of magnetite. For all three surfaces the SX data yield the same 
spectral line shape with a broad maximum at « —0.6 eV and a strong suppression of intensity towards 
the chemical potential. We hence conclude that this is the intrinsic low-energy spectrum of magnetite. 

For a further analysis of the d 6 — > d 5 signal we now turn to its temperature dependence. Figure 
|2Ja) shows SX spectra of the fractured surface measured for temperatures within a few Kelvin around 
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FIG. 2: (a) Spectra of the fractured sample taken around the Verwey transition. Inset: spectral onset energy 
during cooling (down-triangles) and heating (up-triangles) . Also shown is the conductivity hysteresis (solid curves) 
measured on the same sample, (b) Temperature-dependence of the spectral onset energy (see inset for definition) 
over a wider temperature range (see inset for definition). Also shown are the results for the annealed (111) and 
(100) surfaces. 



Ty. Below the Verwey transition the spectrum is consistent with a small energy gap, as expected 
for an insulator. Going through the Verwey transition the spectral onset becomes abruptly shifted 
towards the chemical potential. However, no indication of a metallic Fermi edge or quasiparticle feature 
is observed in the high-temperature phase, in notable contrast to the spectra of oxide materials with a 
Mott transition such as, e.g., V2O3 Hi . For a more quantitative examination we define the spectral 
onset phenomenologically as the intersection of the leading edge with the zero intensity base line (see 
inset of Fig. Hfb)). A plot of this phenomenological parameter versus temperature is shown in the inset 
of Fig. [2Ja). The onset energy jumps exactly at TV and is consistent with the hysteretic behavior of 
the conductivity, confirming that these spectra reflect intrinsic bulk behavior. This conclusion is further 
corroborated by the observation that the onset energy agrees well with the activation energy of the 
conductivity on both sides of the transition • If the onset energy is indeed identified with the insulator 
half-gap (with \x assumed to be in mid-gap position), its discontinuous change AE on » 50 meV will cause 
the concentration of thermally activated charge carriers to increase by a factor of exp(AE on / ksTy) . This 
accounts quantitativel y fo r the observed two-order of magnitude jump in the conductivity, as was already 
argued by Park et al. [lfj . A decrease rather than a complete closing of the gap at Ty also identifies the 
Verwey transition as an insulator-to- insulator transition [lOj . 

We also note that temperature effects are not restricted to the immediate vicinity of Ty- The discon- 
tinuity observed at the Verwey transition is actually superimposed on top of a continuous energy shift of 
the spectral onset over a wide temperature range, as shown in Fig.Efb). The identification of a first-order 
Verwey transition thus requires a very fine temperature grid around Ty. Previous photoemission studies, 
besides being strongly surface-sensitive, used rather wide temperature steps, so that the reported effects 
resulted (at least partly) from the gradual temperature evolution and not from the Verwey transition. 

Interestingly, the discontinuous first-order transition of the onset is only observed for fractured surfaces, 
but not for the (111) and (100) surfaces (also shown in Fig. Etb)). For the (111) sample the onset 
shifts continuously towards the chemical potential and reaches it at ~ Ty without any evidence for a 
discontinuous phase transition or a metallic Fermi edge. The onset of the (100) spectra shows also a 
continuous but much smaller temperature shift. Moreover, the onset remains clearly below the chemical 
potential even up to 190 K, well above Ty. It is tempting to relate this finding to reports of room- 
temperature charge order at the (100) surface 0,0]. These surface effects in the temperature evolution 
can be traced back to the fact that any termination of a bulk magnetite crystal will result in a polar surface, 
which in order to minimize its energy has to undergo a geometric (or electronic) relaxation. The resulting 
atomic reconfiguration may then lead to a suppression of the Verwey transition, similar to the adverse 
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effects of other lattice deformations in the volume induced by defects or hydrostatic pressure 0, 0, H3| • 
For the (111) and (100) samples the thickness of the relaxed surface layer is obviously comparable to 
the relatively large photoemission probing depth and therefore extends over at least one unit cell (lattice 
constant a = 8.4 A). As their preparation involves high temperature annealing, the resulting surface 
configurations are in thermal equilibrium. In contrast, crystal fracture at cryogenic temperatures will 
result in a non-equilibrium surface in which the topmost atomic layer may be disturbed but a full surface 
relaxation is hindered by the lack of sufficient thermal energy. The subsequent layers (accessible by SX 
photoemission) will remain frozen in a bulk-like configuration and thus display the Verwey transition. 

In any case, the above results establish that the broad spectral distribution and the absence of a well- 
defined quasiparticle peak near the chemical potential are intrinsic properties of the electron removal 
spectrum and not a surface artefact. The non-observation of a sharp Fermi edge, particularly above Ty, 
is in clear disagreement with LDA band calculations 0] which predict a metallic state (see Fig. Eta))- 
Furthermore, even if one allows for an insulating state as suggested by the activated conductivity on both 
sides of the transition, the inflexion point of the leading edge of the spectrum (which in a band insulator 
would mark the gap edge) lies at much higher binding energy than expected from the transport (half)gap. 
Rather, the gap coincides with the energy of the spectral onset (as defined above), very reminiscent 
of observations recently made in photoemission spectra of quasi-one-dimensional Peierls systems and 
attributed there to polaronic effects 0, • In this picture strong electron-lattice coupling leads to the 
formation of polaronic quasiparticles, i.e. electron (or hole) excitations heavily dressed by virtual phonons. 
Removing an electron from the coupled system (as in photoemission) results in a spectrum consisting 
of a coherent quasiparticle peak, greatly reduced in spectral weight and renormalizcd in energy, and an 
incoherent background of phonon side bands shifting weight away from the chemical potential. If the 
coupling is sufficiently strong, a pseudogap-like behavior of the spectrum will result. Indeed, for magnetite 
strong coupling to the lattice and the formation of small polarons have already been inferred from other 
properties, e.g., the small carrier mobility |23j . the unusual temperature behavior of the conductivity 
above Ty, and the observation of a mid-infrared polaron peak in the optical conductivity |2 ll l25| . from 
which a polaronic binding energy of e p ~ 300 meV has been derived. 

In order to further substantiate the above qualitative picture we have analyzed the data using the 
theoretical model of Alexandrov and Ranninger j2|| , in which the electron removal spectrum of a system 
with strong electron-phonon coupling can be written as 

oo 2n 

1(e) oc e~ g2 N p (e) + ^ e^ 2 L_N p (e + nuo). (1) 

n=l 

The first term represents the polaronic quasiparticle band (N p (e)) renormalizcd in spectral weight and 
width (i.e. its inverse mass) by a factor e~ 9 , where g 2 — £ p /ujq is a dimensionless electron-phonon 
coupling constant and luq a characteristic phonon energy. The remaining intensity is transferred to 
multiple phonon side bands shifted by nojQ away from /i. For strong coupling (g 2 >> 1) the polaron 
quasiparticle weight is thus exponentially suppressed in the photoemission spectrum, but its binding 
energy still defines the lower edge of the transport gap. According to band theory the (non-coupling) 
density of states (DOS) near the chemical potential consists of two peaks, one right at /i and a second one 
at —0.6 eV (Fig. Ota))- The polaronic effects in the electron removal spectrum have now been modeled 
for the first DOS peak by eq. (1) and for the second by a Gaussian of similar width (2?|]. Figure EJb) 
shows the result for the 180 K spectrum of the fractured sample, with the phonon energy scale loq fixed 
to the highest optical mode (70 meV j2{|) and the theoretical spectrum broadened by ~ 100 meV to 
account for instrumental resolution and a possible residual polaron bandwidth. The best agreement with 
experiment is obtained for a coupling strength of g 2 ~ 5, which accurately reproduces the observed 
pseudogap behavior near the chemical potential and the shift of spectral weight to higher binding energy 
with respect to the band theory DOS. The value for g 2 is consistent with both the polaronic binding 
energy and the effective polaron mass (~ 100 . . . 200m ) derived from infrared spectroscopy |25j. This 
analysis corroborates the qualitative arguments from above for a pronounced electron-phonon coupling 
and formation of small polarons in magnetite. 

The strong coupling to the lattice can be traced back to the fact that Fe 2+ in a high-spin d 6 configuration 
is a Jahn- Teller (JT) ion. In fact, the low-temperature distortions of the oxygen octahedra surrounding the 
formally divalent B-site ions can be understood as JT-active modes . Very recent LSDA+U calculations 
based on the refined structure of Ref. 6] confirm this view. They find that the local occupation of 
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FIG. 3: (a): Near-^t spectrum of a fractured Fe3C>4 crystal at 180 K (i.e. above TV) superimposed on the 
band theory DOS for the high-temperature phase |l5j|. cut-off by a Fermi function for the same temperature, 
(b): Model fit (red curve) of the background-subtracted spectrum (black dots) and its decomposition into the 
polaronic single-particle spectrum of Ref. [2(J and an additional Gaussian (see text for details). "QP" denotes 
the location of the exponentially suppressed polaronic quasiparticle band. 



^g-orbitals on the i?-sublattice follows octahedral distortions as expected from the JT effect, except that 
the actual charge disproportionation is less than ±0.1e, leading to non-integer Fe valencies (2.4+/2.6+). 
As a consequence, the low-temperature phase is characterized by orbital order (OO) and much less so (if 
at all) by CO, consistent with the recent diffraction experiments [1,0 111- It is thus tempting to conclude 
that the physics of the Verwey transition contains elements of a cooperative JT effect, which requires 
additional stabilization by local dd Coulomb interaction [51, |2^| . 

We thus arrive at a picture in which the charge carriers in magnetite are small polarons, presumably 
of the JT type, forming a narrow band of high effective mass. At high temperatures they contribute 
to the charge transport by both coherent tunneling and activated hopping, as originally proposed by 
Ihle and Lorenz [jj. Below Ty the polarons condense into an orbitally-ordered phase (possibly with an 
additional small CO amplitude), whose microscopic origin remains to be revealed. This opens a gap 
in the polaronic quasiparticle band and causes the sudden drop in the conductivity. The polarons are 
much more stable than the OO/CO as inferred from our spectra of the annealed surfaces, for which 
the Verwey transition is suppressed but not the polaronic line shape. Our results demonstrate that the 
Verwey transition in magnetite is not a purely electronic effect and that any microscopic model has to 
account for the important role of electron-phonon coupling in the simultaneous presence of interatomic 
Coulomb interaction. 

We thank A.S. Alexandrov, H. Fehske, S. Fratini, T. Kopp, I. Leonov, A. Loidl, A. Pimenov, and 
D. Vollhardt for helpful discussions. The photoemission experiments were performed with approval of the 
Japan Synchrotron Radiation Research Institute (Proposal Nos. 2002B0104-NSl-np and 2003B0159-NSa- 
np). Funding was provided by the Deutsche Forschungsgemeinschaft through SFB 484. M.S. gratefully 
acknowledges financial support by the Japan Society for the Promotion of Science (JSPS). 



* Electronic address: david.schrupp@physik.uni-augsburg.de 

[1] M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 70, 1039 (1998). 

[2] F. Walz, J. Phys.: Condens. Matter 14, R285 (2002). 

[3] E. Verwey, Nature 144, 327 (1939). 



G 



[4] N. Mott, Metal-Insulator Transitions (Taylor and Francis, London, 1990). 
[5] D. Ihle and B. Lorenz, J. Phys. C: Solid State Phys. 18, L647 (1985). 

[6] J. P. Wright, J. P. Attfield, and P. G. Radaelli, Phys. Rev. Lett. 87 266401 (2001); Phys. Rev. B 66, 214422 
(2002). 

[7] J. Garcia et al, Phys. Rev. B 63, 54110 (2001). 
[8] J. Garcia et al, Phys. Rev. Lett. 93, 156408 (2004). 
[9] A. Chainani et al., Phys. Rev. B 51, 17976 (1995). 
[10] J.-H. Park et al, Phys. Rev. B 55, 12813 (1997). 

[11] Probing depth is defined as three times the inelastic mean free path corresponding to 95% of the photoemission 
signal. 

[12] D. Schrupp et al, Ann. Phys. (Leipzig) 13, 78 (2004). 

[13] Y. Saitoh et al, Rev. Sci. Instrum. 71, 3254 (2000). 

[14] V. I. Anisimov et al, Phys. Rev. B 54, 4387 (1996). 

[15] A. Yanase and N. Hamada, J. Phys. Soc. Jpn. 68, 1607 (1999). 

[16] S.-K. Mo et al, Phys. Rev. Lett. 90, 186403 (2003). 

[17] R. Wiesendanger et al., Science 255, 583 (1992). 

[18] I. V. Shvets et al, Phys. Rev. B 70, 155406 (2004). 

[19] G. K. Rozenberg et al, Phys. Rev. B 53, 6482 (1996). 

[20] J. H. V. J. Brabers, F. Walz, and H. Kronmiiller, J. Phys.: Condens. Matter 11, 3679 (1999). 
[21] L. Perfetti et al, Phys. Rev. Lett. 87, 216404 (2001). 
[22] L. Perfetti et al., Phys. Rev. B 66, 075107 (2002). 

[23] N. Tsuda, K. Nasu, A. Yanase, and K. Siratori, Electronic Conduction in Oxides (Springer, Berlin, 1991). 

[24] S. K. Park, T. Ishikawa, and Y. Tokura, Phys. Rev. B 58, 3717 (1998). 

[25] L. V. Gasparov et al., Phys. Rev. B 62, 7939 (2000). 

[26] A. S. Alexandrov and J. Ranninger, Phys. Rev. B 45, 13109 (1992). 

[27] A Gaussian profile is the limiting shape of eq. (1) for strong coupling. Here we are interested in shape and 
position of the pseudogap behavior near /i which is dominated by the polaron effects in the first DOS peak. 
Therefore a simple Gaussian approximation for the second peak is fully sufficient. 

[28] I. Leonov et al, Phys. Rev. Lett. 93, 146404 (2004). 

[29] H.-T. Jeng et al, Phys. Rev. Lett. 93, 156403 (2004). 



